NBI-HE-96-01 
|Ei^/9601043 
January 1996 
revised February 1996 



Wilsonian Approximated Renormalization Group 
for Matrix and Vector Models in 2<d<4 



Shinsuke Nishigaki I 

The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



Abstract 

Wilson's approximation scheme of RG recursion formula dropping momentum dependence 
of the propagators is applied to large- vector and matrix models in dimensions 2 < d < 4 
by making use of their exact solutions in zero dimension. In spite of apparent dependence 
of critical exponents upon the dilatational parameter p involved by the approximation, 
the exact exponents are reproduced for vector models in the limit p — 0. Application to 
matrix models is then reexamined after the same fashion. It predicts critical exponents 
V = 2/d and rj = 2 — d/2 for the tr$^ matrix model. 



e-mail address: nishigaki@nbi.dk 



1. Introduction 



The study of non-gaussian random matrix models initiated by the pioneering papers 
have been combined with Weingarten's idea of discretized quantum gravity p| and 
yielded a thorough understanding of c < 1 noncritical strings 0. There the 'c = 1 
barrier' [Q, traditionally attributed to the tachyonic nature of ground state of bosonic 
strings, obscures itself among the technical difficulty with the matrix model in dimensions 
d > 1: that it does not allow one to perform angular U{N) integration so that the system 
is no more reduced into free fermionic eigenvalues confined in a potential well. Despite 
the importance of uncovering the nature of phase transition of random surfaces around 
c = 1, several previous attempts proposed to circumvent this difficulty, including the 
Brezin-Zinn- Justin program I, i and the light -cone quantization [0], still fail to provide 
us with a reliable prediction even in the 'planar' large- limit. 

This letter is aimed to present an insight into this long-standing problem, by viewing 
the tr $^ matrix model as a Landau- Ginzburg hamiltonian and exploiting Wilson's treat- 
ment of the renormalization group P). It consists of the following steps to derive an RG 
recursion formula: 

i) to separate into its high/low-frequency parts with respect to an arbitrarily 
chosen mass scale p 

ii) to perform functional integration over the high-frequency field by using an approx- 
imation of propagators -|- r) and —d/ dp'^\p=o — > 1/ (const. + r), as well as trunca- 
tion of induced interactions 

iii) to rescale coordinates and the low-frequency field to pull back the renormalized 
action into the same form as the original action. 

The approximation ii) is equivalent to substituting all the loop integrations by zero- 
dimensional combinatorics, which was readily put in our disposal by ref . . This program 
was previously apphed by Ferretti to a c? = 3 matrix model, although he relied upon 
an assumption of universality for the choice of p which proves incorrect in the sequel. 

In this letter I first apply the program to 0(A^)-symmetric vector models in the large- 
N limit |lT| , which have served as a probe to matrix models due to their resemblance 
to and their simplicity relative to the latter [0, ^ . I shall show an apparent p-dependence 
of the exponents involved by the approximation, and then extract the known exact result 
for 2 < (i < 4* by taking the maximal dilatation limit p — > 0. Next I reexamine the 
application to matrix models and obtain the mass exponent v = 2/d and the anomalous 
dimension rj = 2 — d/2 m the same limit. 

2. Wilsonian approximated RG for vector models 

In this section we consider a Euclidean field theory of an iV-component scalar in 
dimensions 2 < d < A with a cutoff, which we choose as the unit of mass. The action is 
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* Due to the IR divergence this approximation is known to be inapphcable for d < 2. 
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Now we introduce an arbitrary mass scale < p < 1 and separate $ into its low and high 
frequency parts accordingly, 
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We aim to integrate over the high frequency in the large-A^ limit and incorporate its 
effect as a renormalized action of the low frequency p is then to play the role of a new 
cutoff. Substituting $ = $ + into the action (|I|) it reads 

Sm = Sm+a[^,4>] + S[^], 
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The quadratic term separates into S[(^] and S'[0] due to momentum conservation. Inte- 
gration over yields the induced action S'[$], determined by 
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Here (■ ■ ■) denotes an average with respect to the measure Vcpe ^^'^^ 

It is easy to confirm that the last three terms in eq.(^) containing ($ ■ 0) do not 
contribute to S in the large- iV limit. Furthermore we truncate induced interactions to 
those already present in the original action (|ip. This truncation to first three relevant 
terms is a natural approximation for handling the RG transformation, although it is 
justified only a posteriori. Then the terms exhibited explicitly in eq.(^) suffice. By 
reexponentiating the rhs of eq.(^ we obtain 
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(5) 

Now we are in a position to apply the Wilsonian approximation to 0-correlators: to 
replace all propagators l/(p^ + r) appearing in the loop integrals with l/(const. + r). Since 
the final result (eqs.(|l6|,|l8D) is insensitive to the numerical value of the constant, we set 
it equal to unity. This approximation virtually reduces correlators to zero-dimensional 
ones, which are exactly calculable using the saddle point method |[Tl|]. The (contracted) 
two-point function in eq.(|^). 



+ 



■CO + 



-4^ 



d'^p 



N 



d'^q N 
{2txY {p^ + r)2 J {2nY q^ + r 



+ 



(the real/dotted lines stand for 

approx. / Cd 



fields respectively) is approximated to be 
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where q = /p<|p|<i and C2 denotes the zero- dimensional two-point function in 



eg. ( |A6| ) . Similarly the connected four-point function at zero-external momenta is ap- 
proximated using its zero- dimensional counterpart C4 in eq.( |A7| ), 
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A characteristic feature of vector models is that the momentum-dependent self-energy 
diagrams which could have contributed to wave function renormalization vanish in the 
large- A^ limit. Thus the renormalized action for low frequency $ reads 
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with g = u Cd/ (1 + r)^. Note that the coefficient functions of $^ and y^^j by construction 
have an interpretation as normalized IPI vertices [H, 
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l-2gC^{g) 



89 



(9a) 
(9b) 



They can be confirmed by Schwinger- Dyson equations[]; eq.(9B) (xSg) is by 

X. + )(D< = ')&■'. (10) 

(the solid/crossed blobs represent Green/lPI- vertex functions C4/r4 respectively) and 
eq.(|9|) is by eq.(^). 

Finally we must rescale x — >■ p~^x {p ^ pp) and $ — > p'^/'^~^^ so that the renormalized 
action (^ has the same momentum range < p < 1 as the original one (|T|), 

S'm = I d'x [^(VK)^ + ^ + 4(1 + r)g C^] $ ^ + ^p"'' [1 - 2g C,{g)] 
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Therefore the RG recursion equation takes the form [e = 4 — d) 

-1 + (1 + 16^)1/2 
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t I thank G. Ferretti for pointing out this observation. 
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Now we are ready to solve the RG equation following the general scheme. The non- 
gaussian fixed point is determined by m' = u 7^ in eq.(|12b|), that is 



1 + (1 + 16oJ-^/2 
= — ^ — — • (13) 

For e ^ 1 eq. ([T3|) always has a solution ~ e/4 log(l/p) > as it should. For e ~ 1, 
however, it ceases to have a solution on the perturbative sheet (1 + 16 (7)^/^ = + lQg 
for a sufficiently small p because its rhs always exceeds 1/2. Since whether p^2~^/'^ has 
no physical significance we are obliged to continue the fixed point to the second sheet 
(1 + IQgy/'^ = — v^l + 16 g. Consequently the fixed point moves from g = to g = +00 
on the first sheet and then turns back to = on the second sheet as p decreases from 1 
to 0. Under this agreement the fixed point r = r' = r^,, u = u' = u*, g^ = Cd/ (1 + 
in eq.(12) is given by 

9* = ^ 14a 

4(l-2p^)^ 

<0, (14b) 
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p^ (1-p^) (l-p2)' , , 

= — — ^ — ^ > 0. 14c 

4Q(l-p2_pe + 2p2+^)' ^ ^ 

The signs of and u^^ are in accord with the general feature of the Wilson- Fisher fixed 
point for 2 < d < 4. In order to calculate critical exponents we need to linearize the RG 
equation (12) around (r*, m^,), 
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The eigenvalues of the matrix above, 
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Ai,2 = 1 + P''V2 - 3p72 + p2^ ± V(l + p^-72 - 3 pV2 + p^'Y - p2^-2, (16) 
determines the first two of the series of scahng indices i/m. (m = 1, 2, ■ ■ ■) by p 

Vm = — j , (17) 

logp 



the greatest of which is related to the mass exponent 



yi 2 4: \8 2(l-p2); \16 2 (1 - p^ 

Due to the approximation l/(p^ + r) — »• 1/(1 + r), O(e^) or higher order terms depend 
on p, the portion of integrated momentum region by one step and scaling is apparently 
broken. However we can confirm that 7/1,2 are smooth under the switchover of the sheets 
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at p = 2 Moreover, in the p — > limitQthe eigenvalues approach Ai — > ^, A2 — ^ p^ 
and provide us with the exact values for large- vector/spherical models in 2 < d < 4, 



Urn = d — 2m, z/ = l/{d — 2) |jTO|, |l^. We expect that Um for m > 3 be reproduced by 
relaxing the truncation of induced interactions. This exactness might be attributed to 
that in the limit p — the cutoff theory is so strongly course-grained by a single step of 
RG transformation that it flows into the limiting IR theory quickly enough to exceed the 
accumulation of errors in the approximation. We will exploit this observation to calculate 
critical exponents of matrix models in the subsequent section. 

3. Wilsonian approximated RG for matrix models 

Application of the Wilsonian approximation to matrix models was already considered 
in ref. 0; here the outline of derivation of the RG equation is briefly recalled. 

We start from a Euclidean action of an x hermitian matrix field 



S[<^= j d'^xti 



2^ ' 2 N . 
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equipped with a cutoff =1 as before. We separate $ with respect to a momentum p into 
low/high-frequency parts $ = $ + 0, whose coupling reads 



(t[$,</)] = ^ / d^'xti 
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In the case of matrix models, ^-integration also induces products of traces such as ( tr $^)^ 
which is as relevant as the single trace, tr$^. However we can still truncate induced 
interactions to those present in the action (|1^) consistently because the (tr$^)^ term is 
induced always with a suppression factor 1/iV^ relative to tr$^ and thus negligible in the 
iV — s> 00 limit. Taking into account planarity of large-A^ matrix models, Z2 symmetry 
(J) —(j) and momentum conservation, the induced action reads 




m = 

IM J ' 

tr $0^(x) tr $0^(1/)) + ( tr $^0^(x) tr $^0^(i/) \ 

tr <IV^(a;) tr ^<f{y) tr <l0^(z)^ 
tr$0^(x)tr$0^(?/)tr<l0^(^)trl>0^(w)) . (21) 

Again we replace all propagators l/(p^+r) by 1/(1 +r) to approximate mass- and coupling 
constant renormalization. The first two terms in eq. (|2ll) contribute to the tr $^ term. By 
making use of the SD equation, diagrammatically written as (the solid/crossed blobs 
represent Green/ IPI- vertex functions C2/r2,4, respectively) 

Although p ^ corresponds to 0, it should not be confused with that the perturbative 

calculation suffices; the fixed point is indeed on the second sheet of mapping g 1-^ C2,i{g) and this 
double-sheeted structure of Green functions is highly nonperturbative. 
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their contribution are summarized into {g = ucd/{l + r)^) 



5,[$] = (l + r) 



r2(^?) 



(23) 



Similarly, the induced interaction terms (the last three in eq.(^)) are neatly compiled 
into 

'"rf°'xtr$^ 



N \ Ag 



(24) 



under our approximation. 



In the case of large- matrix models, the second term in the Ihs of eq.(22) contributes 
also to wave function renormalization. To incorporate its contribution we need to differ- 
entiate it by the external momentum at p = 0, which can not be treated in the original 
ultra-local approximation. Following Golner's modification |]14| justified on the dimen- 
sional ground, we approximate this procedure simply by replacing with multiplication of 
a propagator 



d 



approx. 



p=0 
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Then the induced kinetic term is evaluated as 



-\ 2 



= 2g V,{g)C2{gY / d'xti V$ 



To recapitulate, the renormalized action for low-frequency $ reads 
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After rescaling the kinetic term to a standard form and then x ^ p ^x, $ 
the RG recursion equation takes the form 
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We are now ready to solve the RG eq. utilizing zero- dimensional Green functions 
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with 12g + a? - 
determined by eq. 



1 = 0. For < e < 2 there exists a unique non-gaussian fixed point 



P 



(30) 



[1 + %, T,{g,)C2{g.ff' 

which turns out to move from ^f^, = to oo as p = 1 ^ 0, always on the perturbative 
sheet a? = (—1 + y/1 + 48(y')/(24(yf)|. We can again confirm < and > 0. The yi^2 
indices are obtained by following the same procedure as in the previous section. For any 
p they turn out to lie in the range yi > > y2 (and are equal to the mean field values 
yi = 2, 1/2 = for e — *• as should be), in accordance with the fact that the Wilson- 
Fisher fixed point for 2 < d < 4 is associated with a single relevant perturbation. The 
O (e^) and higher order terms of the u exponent depend on p as before, and converges to 
z/ ^ 1/(2 - e/2) = 2/d in the p ^ limit (Fig.l). 




Fig.l: Plot of l/i/ for 2 < d < 4. 
Prom top to bottom (at d = 2): p = I, 1/2, 1/4, 10^\ 10"^ 10"^, 10"^ and 0. 

On the other hand, the anomalous dimension rj, determined by the wave function 
renormalization factor via 



_ log [1 + 4(7, T4{g*)C2ig*f] ^ e _ log(r4((?.)/4(7,) 
^ logp 2 21ogp 

can be shown to converge toe/2 = 2 — (i/2in the p — limit (Fig. 2). 



(31) 



^ The wave function renormalization factor is responsible for this fact; without it would proceed to 
the second sheet as in the case of vector models and fail to possess a meaningful p — > limit for critical 
exponents this time. 
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Fig.2: Plot of for 2 < d < 4. 

From bottom to top: p = 1 (r/ = 0), 1/2, 1/4, 10"\ 10"^^ 



10 



-3 



10-"^ and 0. 



Although the exactness of the p — > limit in the case of vector models does not necessarily 
imply that in matrix models, it nevertheless provides us with a strong supporting ground 
for these limiting values of the critical exponents. 

4. Concluding remarks 

In this letter I have identified the Wilson-Fisher fixed point (in the stable m > region) 
both for large-A^ vector and matrix models using the Wilson's scheme and computed 
critical exponents. The essential approximations employed are three-term truncation of 
induced interactions and zero-dimensionalization of the propagators, combined with the 
p ^ limit. This program is, in a sense, complementary to the standard approach []13| 
using the gap equation where the approximation l/(p^ + r) 1/(1 + r) is known to yield 
the full series of exact ym indices for large- iV vector models without ambiguity in p, after 
taking all the induced interactions into account. What is remarkable for matrix models is 
that the critical exponents for magnetization and specific heat, derived from (non- mean- 
field) z/ and 7] via (hyper-)scaling relations, are predicted to stay at the classical mean 
field values a = 0, /5 = 1/2, 7 = 1, 6 = 3, despite d is below the upper critical dimension 
4. This consequence is nontrivial and may not be attributed to the roughness of the 
approximation when we recall its exactness for vector models. 

Generalization to higher-order truncation and criticality as well as to non-hermitian 
matrix models is straightforward. Direct calculation of various magnetic exponents will 
be made possible by relaxing the Z2 symmetry ^ —0, and serve for the check of 
consistency. I hope to discuss these points in a subsequent publication. Application 
of our program to large-A^ QCD utilizing its low-dimensional exact solution is another 
interesting subject, although a special care is required for a cutoff procedure in order to 
maintain gauge invar iance. 

Finally I list a few points yet to be clarified. The precise mechanism for the three- 
term truncation in the p ^ limit to work out successfully for vector models must be 
fully explained in order to justify the matrix model results. Turning back to the original 



8 



motivation, the relationship of these field-theoretic (spacetime) exponents to geometrical 
(world-sheet) exponents d^, etc., measured numerically for c > 1 candidates |T5[, is 
also unclear to the author at present. 
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Appendix: Green functions of OD vector model 

Here I summarize the derivation of Green functions of zero-dimensional vector model 



following ref. The partition function is 



d 



N . 



r dt e-^[^*+5 ^"-5*1= fdt e- 
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where e* 



b'^/N. The saddle point ts in the above is determined by 

1 



The four-point function is 



On the other hand, making use of the SD equation 

^ d 

the two-point function is given by 



0(1). 
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To recapitulate, 
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(A4) 



(A5) 



(A6) 



(A7) 
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